1 


N91-22973 


The Thermal Conditions of Venus 


Vladimir N. Zharkov and V.S. Solomatov 
Schmidt Institute of the Physics of the Earth 


ABSTRACT 

This paper examines models of Venus’ thermal evolution. The models 
include the core which is capable of solidifying when the core’s temperature 
drops below the liquidus curve, the mantle which is proposed as divided into 
two independent of the convecting layers (upper and lower mantle), and the 
cold crust which maintains a temperature on the surface of the convective 
mantle close to 1200°C. The models are based on the approximation ol 
parametrized convection, modified here to account for new investigations 
of convection in a medium with complex rheology. 

Venus’ thermal evolution, examined from the point when gravitational 
differentiation of the planet was completed (4.6 billion years ago), is divided 
into three periods: (1) adaptation of the upper mantle to the thermal 
regime of the lower mantle: approximately 0.5 billion years; (2) entry of 
the entire mantle into the asymptotic regime approximately three to four 
billion years; (3) asymptotic regime. The parameters of a convective planet 
in an asymptotic regime are not dependent on the initial conditions (the 
planet “forgets” its initial state) and are found analytically. The thermal 
flux in the current epoch is ~50 ergs cm’s' 1 . We consider the connection 
between the thermal regime of Venus’ core and its lack of magnetic field. 
After comparing the current thermal state of thermal models of Venus and 
Earth, together with the latest research on melting of the Fe-FeS system 
and the phase diagram of iron, we propose that Venus’ lack of its own 
magnetic field is related to the fact that Venus’ core does not solidify in the 
contemporary epoch. The particular situation of the iron triple point (7 
€ - melt) strengthens this conclusion. We discuss the thermal regime of the 

174 

C-3 


AMERICAN AND SOVIET RESEARCH 


175 


Venusian crust. We demonstrate that convection in the lower portion of 
the crust plays a minor role in regions with a particular crust composition, 
but that effusive or intrusive heat transport by melt, formed from melting 
of the crust’s lower horizons, is the dominant mechanism for heat transport 
to the surface. 

MODIFIED APPROXIMATION OF PARAMETRIZED CONVECTION 

Models of Venus’ thermal evolution, calculated in approximation of 
parameterized convection (APC), were examined in the works of Schubert 
(1979); Thrcotte et al. (1979); Stevenson et al . (1983); Solomatov et al. 
(1986); and Solomatov et al (1987). In parameterizing, dependencies were 
used that were obtained from studying convection in a liquid with constant 
viscosity. They were inferred to be true in cases of more complex rheol- 
ogy. New numerical investigations of convection in media with rheology 
that is more appropriate for the mantle (Christensen 1984a, b, 1985a, b) 
and theoretical research (Solomatov and Zharkov 1989) necessitated the 
construction of a modified APC (MAPC). 

Let the law of viscosity be expressed as (Zharkov 1983): 

tj = 6/r m - 1 exp [A 0 /T ( p/p „) L ] (1) 

where r denotes the second invariant of the tensor of tangential stresses, 
p is density; b and L are considered here to be the constants within the 
upper and lower mantles; m « 3; A,, denotes the enthalpy of activation for 
self diffusion (in K); and where p = p 0 is the reference value of density 
selected at the surface of each layer. 

Tb describe convection in such a medium, following Christensen (1985 
a), we will use two Raleigh figures: 


n agpATd 3 

rtdo — j 

XT)o „ 

(2) 

„ agpATd 3 

Rar = , 

xnr 

(3) 


where a denotes the thermal expansion coefficient; g is the acceleration of 
gravity; AT is the mean superadiabatic temperature difference in the layer; 
d denotes layer thickness, \ is the coefficient of thermal diffusivity; and rj 0 
and rfr are defined by the formulae: 



( 5 ) 
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The dependence of the Nusselt number (determined in terms oHhe 
thickness of the thermal boundary layer 6 ) on Ra<, and Ra Tl where m - 3, 


is parameterized by the formula 


Nu=4i = 

2 b 

(6) 

In the case of free boundaries (lower mantle): 


a = 0.29, fa = 0-37, ftp = 0.16. 

(7) 

In the case of fixed boundaries (upper mantle): 


a = 0.13, Po = 0.35, Pt = 0.15. 

(8) 


The theoretical expression has the following appearance (Solomatov 
and Zharkov 1989): 


Nu = , _ J^tL Rgf^Raf^ = 0.2ZRa° o A Ra° T 2 . (9) 

It is obtained from the balance of the capacities of viscous dissipation and 
buoyancy forces, and is in good agreement with (7) and (8). 

The heat flow at the upper or lower boundary of the layer is equal to 


Fi = 


aeATj 

6 


( 10 ) 


where AT, denotes the temperature difference across the thermal boundary 
layer. Velocities at the boundary (u), mean tangential stresses in the layer 
(r) and mean viscosity (rj ) are estimated by the formulae (m=3). 


7T Xd2 

U — — -r - , 

(ii) 

4 6 

[nub A fp\ L V /3 

(^- eXP T^UJ j ’ 

(12) 

&2 A 0 (p\ L 

^=¥ exp YUj • 

(13) 


The mean temperature of the layer, T , and the temperature of the top of 
the lower thermal boundary layer, T t , can be calculated from the adiabatic 
relationship through the temperature in the base of the upper thermal 
boundary layer, T u : 


T = nT u ; T L = n L T u , 


(14) 


where n and are constants. 
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With large Ra^/Ra^, the very viscous upper thermal boundary layer 
becomes reduced in mobility, taking an increasingly less effective part in 
convection, and the previous formulae are not applicable. Solomatov and 
Zharkov (1989) estimated that the transition to a new convective regime 
occurs when 


Rclt ~ R&Ttr ~ 2R a l] m = 3. (15) 


DESCRIPTION OF THE MODEL 

There is no unequivocal answer to the question of whether convec- 
tion in the Earth’s mantle (or Venus’) is single layered or double-layered. 
Previous works explored the single-layered models of convection. Possi- 
ble differences in the single- and double -layered model of Venus’ thermal 
evolution were discussed by Solomatov et al . (1986) and Solomatov et al . 
(1987). We propose here that convection is double-layered, and the bound- 
ary division coincides with the boundary of the second phase transition at 
a depth of approximately 756 kilometers (Zharkov 1983). 

The thermal model of Venus (Figure la) contains a cold crust, whose 
role is to maintain the temperature in its base at approximately 1200° C 
(melting temperature of basalts), the convective upper mantle, the con- 
vective lower mantle, and the core. An averaged, spherically symmetric 
distribution, T (r,t) is completely determined by the temperatures indicated 
in Figure lb. 

The thermal balance equations for the upper mantle, the lower mantle, 
and the core are written as: 

4 (IT 

-*(Rl - R 3 n )piC pl ^ = 47 rR 2 2 F 12 - 4 xR 2 l F l , (16) 

\*(R?n-&c)p2C P ^ = ^(R 3 2 -R 3 c )p 2 Q 2 -4nR 2 7 F 21 +4*R 2 c F e> (17) 

TnR 3 C pcPc ^ e + Q c ^ = 4 *R 2 c F c . (18) 

Indices “1,” “2,” and “C” relate, respectively, to the upper mantle, 
lower mantle and the core. T denotes the mean layer temperature, Fl is 
the heat flow from the mantle under the lithosphere. The heat flow at the 
surface of the planet is obtained by adding to F^, thermal flow generated 
by the radiogenic production of heat in the crust (~11 erg cm -2 s -1 ). The 
radius of the lithosphere boundary (Rl) differs little from the radius of 
the planet (R„) so that Rj, « R 0 . It is supposed that almost all of the 
radioactive elements of the upper mantle migrated into the crust when the 
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FIGURE 1 (a) Diagram of Venus’s internal structure. Radii are indicated for. the 

planet — R 0 ; the base of the lithosphere-Rz.; the boundary between the upper and 
lower mantles — R 12 ; the core-R c ; and the solid internal core -Ri- The dashed line 
illustrates the boundary of the lithosphere, which is an isothermic surface of T B 
1200° C = const; d cr and denote the thickness of the crust and the lithosphere, (b) 
Schematic, spherically symmetric temperature distribution in the cores of Venus. Reference 
temperatures are indicated for. the surface-T 3 ; lithosphere base-Tj,; base of the upper 
thermal boundary layer of the upper mantle— T ul ; peak of the lower boundary layer of 
the upper mantle— T^i; boundary between the upper and lower mantIes-Ti 2 ; base ot 
the upper boundary layer of the lower mantle— T u2 ; peak of the lower boundary layer 
of the lower mantlc-T^i and the boundary between the core and the mantle— T em . 
Thicknesses of the thermal boundary layers are given: for the boundaries of the upper 

mantle- <5 2 and S c for the boundaries of the lower mantle. The dashed line indicates 
the core melting curve, which intersects the adiabatic temperature curve at the boundary 
between the outer, liquid and inner, solid core. 


crust was melted. F c denotes the heat flow from the core. The heat flow 
at the boundary between the upper and lower mantles with a radius of R l2 
is continuous: F 12 = F 2 i. 

Heat production in the lower mantle (Q 2 ) is defined by the sum 

4 

q 2 = Yj Qoi<* p [M<o - 01 ( 19 ) 

1 = 1 

where Q oi and A, denote the current heat production of the radioactive 
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isotopes K, U, and Th for one kilogram of undifferentiated silicate reservoir 
of the mantle and their decay constants. The concentrations of K, U, and 
Th are selected in accordance with O’Nions el al. (1979): U = 20 mg/t, 
K/U = 10 4 , Th/U = 4. 

The term, Q c dm/dt, in (18) describes heat release occurring when 
the core solidifies after the core adiabat drops below the core liquidus 
curve. It is supposed that the core consists of the mixture, Fe-FeS, and 
as solidification begins from the center of the planet, sulfur remains in 
the liquid layer, reducing the solidification temperature. The value Q c is 
composed of the heat of the phase transition and gravitational energy. 

According to the estimates of Loper (1978); Stevenson el al. (1983), 
and Solomatov and Zharkov (1989), Q c = (1-2) - 10 10 erg g -1 . 

Mean temperature of the adiabatic core: 

T c * n c T CM (20) 

where n c « 1.2, the mean core density is p - 10.5 g cm -3 , C P c = 4.7- 10 6 
erg g -1 K -1 (Zharkov and Tfubitsyn 1980; Zharkov 1983). 

The formulae for the melting T m (p) and adiabatic T ad (p) curves are 
written as follows: 


( x 2.24 


Tad(p) = T cr 


(- 

\Pcm j 


P 

Pcm 


- 1.224 - 0.009405 -^-- 0.1 586 ^ j -0.05672 ^0 . 


( 21 ) 

( 22 ) 

(23) 


Here T 0 denotes the temperature at which pure iron melts at the boundaiy 
core of the radius R c , where p = pcm = 9.59 g cm -3 ; a ss 2; x is the mass 
portion of sulfur in the liquid core, depending upon the radius of the solid 
core R,, and the overall amount of sulfur in the core, x„: 


x(Ri) = 


XoR-c 

R 3 c-Rf 


(24) 


The intersection of (21) and (22) defines the radius of the solidified core 
R, (Figure lb). 

Convection in the upper mantle is parameterized by the MAPC with 
the parameters (7) (for fixed boundaries) in the lower mantle; and by the 
same formulae with the parameters (8) (for free boundaries). However, 
the difference between (7) and (8) is not very substantial. The parameters 
for the upper mantle are: bj = 4.3 ■ 10 15 dyne 3 cm -6 s, A«i = 6.9 • lff'K, 
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ocj = 3 - 10- 5 K-\ p x = 3.7 g cm- 3 , X i = lO-’cnrV 1 , ®, = 4.5 • 10 5 erg 
cm _1 s- 10 7 erg g 

and for the lower mantle: 

b 2 = 1.1 • 10 17 dyne 3 cm- 6 s, = 1.3 • K^K, <x 2 = 1.5 • lO-'K" 1 , p 2 

= 4.9 g cm" 3 , X2 = 3 • 10- 2 cm _2 s _1 , ae 2 = 1.8 • 10 6 erg cm^s ! K \ g = 
900cm s~ 2 , Cp 2 = 1.2 • 10 7 erg g" 1 , n 2 = 1.13, n L2 = 1-26. 

NUMERICAL RESULTS AND ASYMPOTOTIC SOLUTION OF THE 

MAPC EQUATIONS 

The evolution of the planet began approximately 4.6 billion years ago. 
The initial state of the planet was a variable parameter. The initial value 
of T u2 is the most significant, since T„i, due to the low thermal inertia 
of the upper mantle, rapidly adapts to the thermal regime of the lower 
mantle (t < 0.5 billion years). The core has little effect on the evolution of 
the planet in general and the majority of models do not take into account 
its influence. T u2 was selected as equal to 2500, 3000, 3500 K. The upper 
value is limited by the melting temperature of the mantle, since a melted 
mantle is rapidly freed from excess heat. 

The upper mantle adapts to the thermal regime of the lower mantle for 
the first approximately 0.5 billion years (Figure 2a, b). Then, after approx- 
imately three to four billion years the entire mantle enters an asymptotic 
regime which is not dependent upon the initial conditions. The evolution 
picture in general is similar to the one described by Solomatov et al. (1986); 
and Solomatov et al. (1987). The planet is close to an asymptotic state in 
the present epoch. Contemporary parameters of the models are: 

T ul = (1700-1720) K, U! = (2.1-2.4)cm yr" 1 

Ti 2 = (2500-2530) K, u 2 = (0.8-1.0)cm yr 1 

T u2 = (2840-2870) K, n = (5-6) bars, 

Fi = (35-40)erg cm" 2 s" 1 , r 2 = (110-120) bars, 

<n = (28-30)km, m = (I^IO 21 poise, 

a 2 = (130-140)km, = (3-10)10 22 poise. 

According to the criterion (15) MAPC are applicable throughout the 
entire evolution, just as with the quasistationary criterion (Solomatov et al. 
1987). 

The asymptotic expression for F L in the first approximation is formu- 
lated as (Solomatov et al. 1987): 

f l = ( 1 + ) . (“) 

where Fq denotes the thermal flow created by the radioactivity of the lower 
mantle, and 4 is the characteristic time of decay: 
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FIGURE 2 (a) Evolution of the base temperatures, T u i, T\i and T U 2> with differing 
initial conditions. The core is not taken into account and is considered to be T cm * T/^2 
(Figure 1). (b) Evolution of the thermal flow under the lithosphere, F/,, with differing 
initial conditions. The thermal flow to the surface is obtained by adding ~ 11 erg cm^s” 1 
to Fl* generated by the radioactive elements of the crust. The dashed line indicates the 
thermal flow generated ty radioactive elements of the lower mantle. 
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equal ~5.67 • 10 9 years at the current time. 

The time scale of thermal inertia of the mantle is equal to 

tin = tinl + <ir»2(l + ^)> 

where t, nl and t, n2 are the time scales for the inertia of each of the layers 
individually: 

l,„, = (l + ft + * A ' ( % - ^ )" , (28) 

tin7 = (i+/? 2 + T ^) , (29) 

ATi 1 + l)).|8 1 + 00 + 02 - 1 /3(/?o - W7)Aq 2(T u 2 - Tu)T~2 (jq-j 
^ = AT^ 7 l + Pi+fiMTui-TBYIZi 

Mi and M 2 denote the layer masses, s = R 12 2 / Rl 2 - The values of (27) 
through (30) are calculated in a zero approximation: F L = F Q- and A = 

A<,/n ( p/p 0 ) L • J . ^ ... . 

The time scale for thermal inertia, as compared with models based on 
the conventional APC (Solomatov el at. 1986, 1987), has increased from 
~ 2.5 • 10 9 years to ~3.5 • 10 9 years. The mantle and core temperatures 
obtained are somewhat lower (by 300 K near the core). In the models with 
the core, T em » 3720K and F c w 15 erg citrV 1 in the contemporary 
epoch. Since the adiabatic value is F c « 30 erg cm 2 s there is no con- 
vection if solidification is absent, and a magnetic field cannot be generated 
(Stevenson et al. 1983). 

MAGNETISM AND THE THERMAL REGIME OF 
THE CORES OF EARTH AND VENUS 

We will estimate the temperatures in the Earth using the MAPC. 
MAPC is not applicable to the Earth’s upper mantle, since convection in 
the Earth involves the surface layer, and rheology is, in general, more 
complex. However, MAPC is applicable to the Earth’s lower mantle. Let 
us assume a value for the temperature at the boundary between Earth s 
upper and lower mantles of 

Tj2 = (2300 - 2500) K, (31) 

(Zharkov 1983), and the thermal flow at this boundary is 
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F 2 \ « Jjr (f - F cr + CpMi^^j « 80 ergcm 2 s l , 

where F - F cr w 70 erg cm -2 s -1 is the medium thermal flow from the 
Earth, after deducting heat release in the crust (according to the estimates 
of Sclater et al. 1980), dTi/dt sa -lOOK/billion years. (Basaltic Volcanism 
Study Project 1981). 

We will assume the following values for the rheologic parameters: 


b = 6.6 • 10 16 dj/ne 3 cm _6 s, A = 1.5 • 10 5 K, andA 0 = 1.3 • 10 5 K. 

The remaining parameters are listed in Zharkov (1983). 

As a result, we have: 

T u2 = (2800-2900) K, 

T L2 = (3600-3800) K, 

T cm = (3800-4000) K, 

S 2 = 100 km. 

The value of T cm - T L2 depends upon F c , which we will assume to be equal 
to the adiabatic value of F a< t « 30 erg cm~ 2 s _l . 

Therefore, the temperature at the boundary of the Earth’s core, T C me 
is greater than for Venus (T cmv) by a value of 

Tcme — Tcmv = (100 — 300)/f. (33) 

In order for Venus’ core not to solidify, it is necessary that the adiabat of 
the Venusian core not drop below the solidification curve during cooling. 
Otherwise, solidification of the core will cause the core to mix by chemical 
or thermal convection, and it will trigger the generation of a magnetic field 
(Stevenson et al. 1983). We will estimate the difference between T cme and 
the temperature, T crV (which is critical for the beginning of solidification), 
at the boundary of Venus’ core, with a single pure iron melting curve, T m 
(P) and a single equation of the state for iron p( P). 

T cme is found from the intersection of the adiabat of the Earth’s 
core (22) and the liquidus curve (23), with a sulfur content in the core 
of xe at the boundary of the Earth’s inner core. T cr v is obtained from 
the intersection of the adiabat of Venus’ core (22) and the liquidus curve 
(23) with a sulfur content of x^. We then obtain (Solomatov and Zharkov 
1989): 


Tcme - T cr v = (+300) + (-300), K, 


(34) 
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where x E - x v = 0 4- 0.07. Therefore, if x v £ x E - (0 - 0.02) = 0.07 
^ 0.12; (with x E = 0.09 4 - 0.12; Aherns 1979), the core of Venus is not 
solidifying at the present time, and a magnetic field is not being generated. 

Complete solidification of the core would have led to the absence 
of a liquid layer in the core and would have made it impossible for the 
magnetic field to be generated. However, for this, the temperature near 
the boundary of Venus’ core should have dropped below the eutectic value, 
which, according to the estimates of Anderson et al. (1987) is ~ 3000 
K and according to Usselman’s estimate (1975) is ~ 2000 K. Such low 
temperatures of the core seem to be of little probability. 

Pressure in the iron triple point, 7 — c— 1, (liquid) approaches the 
pressure in the center of Venus. This is significant for interpretation the 
absence of the planet’s own magnetic field. According to the estimates 
of Anderson (1986), P, p « 2.8 Mbar (Figure 3), but is also impossible 
to rule out the larger values of P (p . In the center of Venus, P cv ~ 2.9 
Mbar; at the boundary of the solid inner Earth’s core, P i E = 3.3.Mbar; 
and in the Earth’s center, P C e = 3.6 Mbar (Zharkov 1983). If P ct) £ 
P tp , the conclusion that there is no solidification of Venus’ core is further 
supported, since, in this case, the core’s adiabat (critical for the beginning 
of solidification) drops several hundred degrees lower. This stems from the 
fact that reduction in the temperature of solidification of the mixture, 7 - 
Fe - FeS is greater than for the mixture c- Fe — FeS by ASJAS^ times, 
where AS f /AS 7 denotes the ratio of enthropy jumps (during melting) equal 
to ~ 2, according to Anderson’s estimates (1986). 

The latest experimental data on the melting of iron allow us to estimate 
the melting temperature in the cores of Earth and Venus. Figure 3 shows 
melting curves obtained by various researchers, and the 7 - e boundary, 
computed by Anderson (1986). With x = 0.09-0.12 and a = 1 4 - 2 (formula 
21), the full spread of temperatures at the boundary of the Earth’s core, 
leading to the intersection of the liquidus and adiabat curves, is equal to 

T C me = 3500 4 - 4700 Jf, (35) 

- with data from Brown and McQueen (1986) and 

Tcm e = 4300 4 - 5400/C, (36) 

- with data from Williams et al. (1987). 

The effect of pressure on viscosity of the lower mantle reduces the 
effective Nusselt number and increases the temperature of T C me and 
Tcmv by ~ 300 K (Solomatov and Zharkov 1989). We obtain the estimate 

Tome = 3800 4 - 4300 K, (37) 

which is the best fit with (35). 
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Pressure, Mbar 

FIGURE 3 Pure iron melting curves according to various data and the boundary of the 
phase transition 7 — Fc — e — Fe. The figures indicate melting curves: 1 is from the study 
by Williams a at. (1987); 2 is from the study by Anderson (1986) and experimental data 
from Brown and McQueen (1986); and 3 is based on the Lindeman formula with the 
Grunehaiscn parameter of G = 1.45 = const, using experimental data from Brown and 
McQueen (1986). The transition boundary, 7 — e, was constructed in the study by Anderson 
(1986) and together with the melting curve gives us the location of the triple point. The 
vertical segments illustrate errors in determining temperature for curves 1 and 2. The 
boundary labels are as follows: CM is the boundary between the core and the mantle; C 
denotes the center of the planet; I is the boundary of the solid inner core; and the final 
letter indicates Earth (E) or Venus (V). 7 - € - 1 denotes the location of the triple point. 


TIIE THERMAL REGIME OF THE VENUSIAN CRUST 

What arc the mechanisms by which heat is removed from Venus’s 
interior to the surface? Global plate tectonics, like Earth’s, are absent on 
the planet (Zharkov 1983; Solomon and Head 1982), although in individual, 
smaller regions, it is possible that there are features of plate tectonics (Head 
and Crumpler 1987). A certain portion of heat may be removed by the 
mechanism of hot spots (Solomon and Head 1982; Morgan and Phillips 
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1983V Conductive transport of heat through the crust clearly plays a large 
role However, with a large crust thickness ( > 40 km), as indicated by 
data from a number of works (Zharkov 1983; Anderson 1980; Solomatov et 
al. 1987), the heat flow of F w (40-50) erg cm 2 s 1 triggers the melting of 
the lower crust layers and removal of approximately one half of heat flow 
by the melted matter. The latter flows to the surface as lava or congeals m 

the crust as intrusions. . 

We shall discuss the possibility of solid-state convection in the crust as 

an alternative mechanism of heat removal through the Venusian crust. 

The crust is constituted of an upper, resilient layer with a thickness of 

d e and a viscous one with a thickness of d„ = d cr - d e . 

The crust thickness of d cr is constrained in our model by the phase 
transition of gabbro-eclogite, since eclogite, with a density higher than 
that of the underlying mantle, will sink into the mantle (Anderson 1980, 
Sobolev and Babeiko 1988). The depth of this boundary is dependent 
upon the composition of basalts of Venus’ crust. We assume that d er - 
70 km (Yoder 1976; Zharkov 1983; Sobolev and Babeiko 1988). It is 
considered that radioactive elements are concentrated in the upper portion 
of the crust, and the primary heating occurs via the heat flow from the 
mantle of Fl « 40 erg cm -2 s -1 . The boundary (d e ) of the resilient crust 
is defined as the surface of the division between the region effectively 
participating in convection and the nonmobile upper layer. Crust rheology 
is described by law (1) with the parameters from Kirby and Kronenberg 
(1987). Four modeled rocks are considered: quartz diorite, anorthosite, 
diabase and albite. Parameter m in (1) is close to 3 for them, so that the 
formulae MAPC with parameters (8) for fixed boundaries is fully suitable 
for estimating. In addition, criterion (15) is used to define the boundary, 
d e . At this boundary, which is regarded as the upper boundary of the 
convective portion of the crust, the temperature is equal to 

T 0 = 733 + 2 Od e (km), K. (38) 

The following physical parameters are assumed (Zharkov et al. 1969). 


— 5 ts — 1 


p = 2.8# cm 3 ,a = 2-10 5 K 


ae 


-i 


= 2 • 10 5 er<jf cm l s 1 K 


Computations have demonstrated that the thickness of the resilient 
crust is 20-30 kilometers. The mean temperature of the convective l^er 
of the crust has been calculated at ~ 1600K, 1700K, 1900K and 2000K, 
respectively, for quartz diorite, anorthosite, diabase, and albite, and exceeds 
the melting temperature for basalts by hundreds of degrees. This means 
that convection does not protect the crust from melting, and heat is removed 
by the melted matter. 
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We can estimate that portion of the heat which is removed by con- 
vection. For this, let us assume the temperature in the base of the crust 
to be T b = 1500K (approximately the melting temperature). We obtain 
the Nusselt number from the formulae of MAPC, and find that it is ~1.7; 
1; 0.6; and 0.5, respectively, for quartz diorite, anorthosite, diabase, and 
albite. For rj = const, and for complex rheology, as in (1), (Christensen 
1984a and 1985a), convection begins at Nu ~ 1.5-2, as calculated according 
to the formulae of APC (MAPC). Thus, where d cr « 70km, perhaps only 
quartz diorite is convective, removing 25-30 erg cm _2 s _1 , while 10-20 erg 
cm _2 s _1 is removed by the melt. The rate of convective currents is three 
to five millimeters per year. It can be demonstrated that 

Nu ~ ( d cr - d e )° 9 (39) 

With d cr ~ 120 km, convection also begins in anorthosite, while in 
quartz diorite it occurs virtually without melting. 

These estimates show that convection in the crust may play some role 
in individual regions, depending on the thickness and composition of the 
crust It is not excluded that in individual regions, convection in the crust 
makes its way to the surface. The bulk of the heat is, apparently, removed 
by melted matter. The rate of circulation of material from the crust is, 
with this kind of volcanism, 50-100 km 3 yr _1 . This is three to five times 
greater than crust generation in the terrestrial spreading zones. Another 
process by which basalt material circulates is where new portions of melted 
basalt reach the crust from the upper mantle, and basalt returns back to the 
mantle in the eclogite phase. This process may trigger the accumulation of 
eclogite in the gravitationally stable region between the upper and lower 
mantles. It may lead to the chemical differentiation of the mantle. This 
process has been noted for the Earth by Ringwood and Irifune (1988). 

Figure 4 illustrates various processes involved in heat and mass trans- 
port. 

The geological structures observed on Venus may be related to these 
processes. Flat regions may be tied to effusive, basalt volcanism. Linear 
structures in the mountainous regions may be related to horizontal defor- 
mations which are an appearance of convections in the crust or mantle. 
Ring structures may stem from melted intrusion or hot plume lifted towards 
the surface from the bottom of the upper mantle. 

Convection both in the mantle and in the crust is, apparently, nonsta- 
tionary (as on Earth). This nonstationary nature stems from instabilities 
occurring in the convective system. The characteristic time scale for such 
fluctuations is t ~ d/u, where d ~ lOMO 9 cm is the characteristic dimension, 
and u ~ 1 cm per year is the characteristic velocity. Therefore, character- 
istic “lifespans” for various occurrences of instability are t ~ 10M0 9 years. 
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FIGURE 4 Process of heat transport through the Venusian crust: (a) transfer of heat 
carried from the mantle by the mechanism of conductive thermal conductivity, (b) heat 
removal by the melted matter which is formed as the basalt crust melts. The melted matter 
may flow to the surface or form intrusions; (c) convection in the crust which does not 
reach the surface or reaching the surface; (d) lifting of hot plume from the bottom of the 
upper mantle to the crust of Venus, triggering enhanced heat flow, a flow of crust material, 
and the sinking of the crust; (e) involvement of the crust in mantle convection with the 
formation of spreading zones; (f) basalt exchange between the crust and the mantle. Basalt 
is formed when the upper mantle partially melts and is returned back as eclogite. 


Regional features of tectonic structures, thermal flows, volcanic activity and 
so on, can exist for this length of time. 


CONCLUSION 

1. Modification of the approximation of parameterized convection 
for the case of nonNewtonian mantle rheology led to no marked increase 
in the time scale of thermal inertia of the mantle from two to three or three 
to four billion years in comparison with the usual parameterization. The 
thermal flow at the surface of Venus of - 50 erg cm“ 2 s _1 , is the product 
of radiogenic heat release from the mantle (50%), heat release in the 
crust (20%), and cooling of the planet (30%). These figures for Earth are 
approximately 40, 20, and 40% for the double-layered convection models. 
Temperatures in the upper portion of the mantle are approximately 1700 K, 
which is 50-100 K greater than on Earth. Given the existing uncertainties 
in the concerning parameters, T cm ~ 3700-4000 K at the core boundary 
and may, possibly, be greater. This temperature is 100-300 K less than for 
the Earth. 

2. The magnetic field on Venus is absent. This is most likely due to 
the lack of core solidification and, respectively, the lack of energy needed 
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to maintain convection in the liquid core. For this, it is enough for the 
sulfur content in Venus’ core to be even a little bit less than on the Earth 
(but not more than 20-30% less). This conclusion is stronger if the triple 
point in the phase diagram for iron, y - e— 1, lies at pressures that are 
greater than in Venus’ center, but approximately less than in Earth’s center. 

3. At the values obtained for heat flows from Venus’ mantle, the crust 
melts, and the melted matter ( < 100 km 3 per year) removes about one half 
of the entire heat. The remainder is removed conductively. In individual 
regions, depending on the crust thickness and composition, a portion of the 
heat may be removed by convection, reducing crust temperature and the 
portion of heat removed by the melted matter. Flow velocities comprise 
several millimeters per year. Due to the insufficiently high temperature of 
the surface, convection in the crust separates from the surface as a highly 
viscous, nonmobile layer with a thickness of 20-30 kilometers. In individual 
regions, crust convection may emerge to the surface. Basalt circulation also 
occurs by another way: the basalt is melted out of the upper mantle and 
returned back in the form of eclogite masses, which sink into the lighter 
mantle rock. It is possible that this process triggers the accumulation of 
eclogite at the boundary between the upper and lower mantles, resulting 
in the chemical separation of the mantle. These processes may explain the 
formation of various geological structures on Venus. 

4. The nonstationary nature of convection in Venus’ mantle and crust 
determine regional features of tectonic, thermal and volcanic appearances 
on the surface of the planet which have a characteristic duration of ~ 
lOMO 9 years. 
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